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O ■ Abstract 

, In heavy-ion collisions the created quark-gluon plasma forms a quickly evolving background, 

leading to a time dependent radiative behavior of high momentum partons traversing the medium. 
We use the Schwinger Keldysh formalism to describe the jet evolution as a non-equilibrium process 
including the Landau-Pomeranschuk-Migdal effect. Concentrating on photon emission, a compari- 
son of our results to a quasistatic calculation shows good agreement, leading to the conclusion that 
the radiative behavior follows the changes in the medium almost instantaneously. 
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I. INTRODUCTION 



Relativistic heavy-ion collision experiments, as performed at the Schwerionen Synchrotron (SIS) 
at the Gesellschaft fiir Schwerionenforschung (GSI), the Super Proton Synchrotron (SPS) at CERN, 
the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) and the 
future Large Hadron Collider (LHC) at CERN, allow for studying strongly interacting matter 
under extreme conditions. One main objective of these experiments is the creation and exploration 
of a new state of matter of deconfined quarks and gluons, the so called quark-gluon plasma (QGP). 
The theory of quantum chromo dynamics (QCD) predicts that it is created at high densities and 
temperatures which occur during heavy-ion collisions at high energy. 

The lifetime of the QGP is of the order of up to 5 fm/c. After that, it transforms into a 
gas of hadrons. Thus experiments cannot directly access the quark-gluon plasma, which makes 
the determination of the properties of this state difficult. Therefore, it is important to find the- 
oretical signatures which provide a distinction between a hadron gas and a quark-gluon plasma. 
Furthermore one has to find experimental observables from which one can draw conclusions on 
these theoretical signatures 0]. An important category of these observables are hard probes, 
i.e., partons with high transverse momenta produced during the early stage of a heavy- ion collision. 
These so called parton jets are significantly modified as they traverse the QGP. The properties of 
the medium can then be inferred by comparing the resulting hadron spectra to those for a simple 
p+p collision where medium effects are absent. Upon comparison to theoretical predictions for a 
hadron gas and a QGP one can draw conclusions on the latter one. The output of RHIC experi- 
ments 0, 13) S] shows that strong jet quenching occurs in heavy-ion collisions. The yields of high 
Pt hadrons are suppressed with respect to the p+p reference by a factor of five. 

Over the past decades, several theoretical formalisms 0,0)S)l9) l3) U) 12] have been established 
to describe energy loss due to bremsstrahlung processes. These radiative processes turned out to 
contribute strongly to the energy loss in a partonic medium [ol. [l^. . while collisional energy 



loss may not be negligible [16|, 113, [la, ll9|, 
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22, 



Analogously to Bremsstrahlung in QED where photons are emitted, a parton which moves 
through the medium is expected to emit gluons. Here it is of particular interest how the standard 
pQFT cross section is modified by multiple scattering of the initial parton 2J, |25|] . In the frame- 
work of QED, this has already been studied by Landau, Pomeranschuk 20] and Migdal 27|, who 
discovered a destructive interference phenomenon of the radiation from different scattering points. 
It occurs at small transverse photon momenta and thus leads to a reduction of the pQFT result in 
that range. The general arguments leading to this so called Landau-Pomeranschuk-Migdal (LPM) 
effect indicate a comparable behavior of radiative processes in QCD. Calculations by Gyulassy and 
Wang 28] as well as subsequent works [g. 
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12l . 1291 ] have indeed verified this. 



Furthermore, the jet traverses a dynamic medium that expands and cools down. Thus the 
density of scattering centers is reduced which makes the radiative behavior of the jet time dependent 
as well. In this work, we study whether it adjusts to changes of the medium instantaneously 
(Markovian case) or if certain memory effects to be described below are of importance. 

The non-equilibrium dynamics of heavy-ion collisions has always been a major motivation for 
investigations on non-equilibrium quantum field theory. In this context, the role of memory effects 
during time evolution has first been investigated by j3d. l3l|. In these works, non Markovian effects 
have been found to appear in the collision term of quantum transport equations. Later on, these 
effects have also been observed in the time evolution of so called disordered chiral condensates [32 1 
as this evolution is governed by dissipative and non Markovian stochastic equations of motion. 
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Furthermore, memory effects also occur in the thermahzation of out of equihbrium (^^-theory 



33, 13J] as well as in the dilepton production from non equilibrated hot hadronic matter j35l . l36l ] . 



In principle, the underlying field theoretical description always leads to equations of motion being 
non-local in time. This indicates that memory effects might also be relevant for the radiative 
energy loss of a jet transversing a dynamic QGP. 

In this work, we use the non-equilibrium Green's function formalism first developed by Keldysh 
and Schwinger 37|, |3a, |39|, |40(] to study this non-equilibrium phenomenon. We approach the problem 
by first considering photon emission. In particular, our work has two main goals. First, we want 
to find the timescales in which the photon production rate adjusts to changes in the medium and 
determine their dependence on the momenta of the emitted photons. Second, we dynamically 
calculate the energy loss of the fermion jet over the time interval in which the plasma expands and 
cools down. Comparing our results to a quasistatic calculation determines the possible importance 
of memory effects for in-medium photon emission. 

This paper is organized as follows. In Section II, we derive an expression which relates the 
photon rate with the current-current-correlator of the emitting system. This expression is non- 
local in time and hence accounts for the finite memory of the system. As source particle, we 
consider an ultrarelativistic fermion jet. In Section III, we present numerical results. After a short 
description of how the time dependencies have been implemented, we determine how the memory 
times of the partial photon rates depend on the momenta of the emitted photons. After that, we 
calculate the memory times for the total photon number and the total radiation power. Finally, 
we consider the energy loss of the fermion jet over the time interval in which the plasma expands 
and cools down. In Section IV, we close with a short summary of our results and an outlook to 
possible future research projects. We also touch on the differences between the QED and QCD 
frameworks. 



II. NON EQUILIBRIUM PHOTON PRODUCTION 



Since we concentrate on the investigation of time dependencies, it is sufficient to consider a 
system which is spatially homogeneous. In this case, the inclusive photon spectrum reads 



A 



(1) 



A 



where N\^H{k,t) denotes the number operator for photons with polarizations A and momentum 
k. Note that operators are indicated by They are expressed in the Heisenberg picture which is 
denoted by the subscript 'H'. As we consider a system which is out of equilibrium, the average is 
taken with respect to the initial density matrix. The creation and annihilation operators can be 
expressed in terms of the electromagnetic field operator and its time derivative: 



ax,H{k,t) 
^x,Hik,t) 



(k) 



/2k 
eW*(fc) 



2k 



iAH{k,t) + kAuik^t) 
iA\j{k,t) - kA''jj{k,t) 



Upon insertion of (|2ap and (|2bp into ([T|), we find 

1 



(2vr)= 



-ff'^ik) {dtdf + ik{dt - df) + k^) iD<^{k, t', t) 



t=t' 



(2a) 
(2b) 

(3) 
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where we have introduced the photon propagator iD^^{k,t' ,t) which is given by 

iD<^ik,t',t) = {Al^{k,t)A,^H{k,t')). 
^^^^{k) denotes the polarization tensor which reads: 



(4) 



A=l 



Here we have introduced the timehke unit vector = (1,0,0,0) as well as the spacelike unit 
vector: 

k^' - (k ■ n)n^' 



^{k-nf-k"^ 



(5) 



Up to here, our calculations have basically followed [41[. In the next step, we have to calculate the 
photon propagator. Instead of employing a strictly perturbative analysis as carried out in [41[, we 
evaluate {k,t' ,t) using a generalized fluctuation-dissipation-relation 



35, 



42, 



D<^ik,t',t) 



t' ft 



dti dt2D^^{k, t', h)U<f^{k, hM)D^^{k, t2, t) 



(6) 
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with the retarded and advanced photon propagators D^/^ and the photon self energy ill^ . Since 
photons are expected to leave the medium undisturbed, D^^^ can be approximated by the free 



propagators D^^^: 



DXoik,t,t') = -- sin {k{t - t')) Q{t - t')gf,, , 
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DXoik,t,t') = - sin (A;(t' - t)) e(t' - % 



(7a) 

(7b) 



Upon insertion of (|7ap -(j7ap into ([3]) and taking the time derivative, we obtain for the photon 
production rate at time t: 



1 



(2vr) 



n^'^ik) 



du ■ 2Re 



(A;,t,M)e 



ik(t—u) 



(8) 



This expression directly relates the photon production rate at a given time t to the full history of the 
emitting system described by the current-current-correlator. This non-locality in time accounts for 
possible memory effects. Equation ([8]) is a generalization of the corresponding equilibrium formula 



4J, |45|], which it reproduces for stationary systems. 

In the case of photon emission from a fermion jet moving through hot and dense matter, 
zn^^(A:, u), the self energy emerging from the fermion-photon interaction, is given by the irre- 
ducible part of the current-current-correlator in coordinate space 



exp 



irrcd 

iytiHi{h)]j^,i{x)ji,{x') 



(9) 



irrcd 



where /*° describes the Keldysh countour integral starting and ending at time to and running 
above the largest time argument of ill^i,{x,x'). Tq implies time ordering of the operators along 
the contour and Hi{ti) denotes the interaction part of the QED-Hamiltonian 



Hiiti 



d^xij^j{xi)A^{xi) . 



(10) 
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The subscript T denotes that the operators are expressed in the interaction picture. In order to 
keep the notation more shorthand, we shall omit this subscript from now on. The current operators 
jf^{x) and jl{x') are given by 

jf,{x) = ei}{x)-i^'4}{x) , (11a) 

jl{x') = ei^{x')^^i,{x'). (lib) 

Irreducibility means that one includes only those diagrams which cannot be separated by cutting 
a single fermion line. Since the expectation value of an odd product of operators always vanishes, 
we can introduce an effective interaction operator 



efr 



d^Xi / d^X2ja{xiW{xi)jp{x2)AP{x2) , 



and rewrite Equation ([9]) as follows: 



n=0 



(2n) 



irrcd 



(12) 



(13) 



Next, we apply a finite temperature Wick decomposition [42t | to each term of (USD . In order to 
avoid double counting, we have to take into account only those diagrams which represent real 
interactions and are not already included in the full fermion propagators. Since the n-th order 
term of (jlSp formally scales with ctg"*"^, it is sufficient for our purpose to consider the n = term 
only, which leads to the standard one-loop-approximation of Eq. ([9|) 

ill^u{x,x') = (Tc j^,ix)jl{x') 

= ec{t,t') {j^,{x)jl{x')) + @c{t',t) {jl{x')j,{x)] 

= ec{t,t')iU>,{x,x') + ec{t\t)iU<,{x,x') , (14) 

Here @c denotes the step function generalized to the time contour C. in^^(a;,x') and il[^^,{x,x') 
are given by 

iU>^{x,x') = (^j^{x)jl{x')^ , (15a) 
iU<,{x,x') = (^jl{x')j^{x)) . (15b) 
Insertion of (|llap - (|llb|) into ()15bp and application of the finite temperature Wick theorem yields: 



ill^^{x,x') = f? • Tr {7^5'^(x,x')7i/'S'^(2;',a;)} . 
Here S'^ix^x'^ and 5^(x',x) are the fermion propagators which are given by 

-^^S^ii^x.x') = {^ji{x')i)a{x)'^ 

iS<p{x',x) = l'4)a{x')i)p{x) 



For a spatially homogeneous system, ill'j2,y{k,t,u) reads 



(A;,t,u) = 



(2^ 



Tr ^-f^S^{p + k,t,u)jyS^{p,u,t)'j . 



(16) 

(17a) 
(17b) 

(18) 



Assumingthat the emitting system is in local equilibrium, the Kubo Martin Schwinger (KMS) 



relation 



371 . kSj l yields for the propagators in frequency space 

S^{p + k,E) = 27rinF{E)AF{p + k,E) , 
S>{p,E) = -27Ti{l-nFiE))AF{p,E) 



(19a) 
(19b) 
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Here Af(p,E) is the spectral function and np[E) denotes the equihbrium Fermi distribution. 
Since we want to consider the radiation arising from a highly energetic fermion jet, we make the 
following ansatz. We associate S^{p + k,E) with the ingoing particle and replace npiE) by an 
effective distribution f{E,p) which includes the fermion jet as one additional particle with fixed 
momentum p': 

npiE) ^ f{E,p + k) = np{E) + ^^S^^\p + k-p') . (20) 

Here V denotes the volumes of the emitting system. Since we are only interested in the radiation 
from the fermion jet, we neglect the Fermi distribution np{E), and S'^{p,E) thus reduces to: 

S<{p + k,E) =2TTi^^6^^\p + k-^)AF{p,E) . (21) 

If we further associate {p, E) with the outgoing particle and assume that the mean occupation 
number of this state is negligibly small before the emission process, we can replace 1 — npiE) by 
unity and obtain 

S>{p,E) = -2mAF{p,E). (22) 

The latter assumption is legitimate as we consider a highly energetic fermion jet scattering into 
unoccupied states. For Af{p,E) we use the fermionic Breit Wigner structure as done in [471]. It 
generalizes the standard onshell energy momentum relation 

Af%p, E) = {Ejo -j-p + m) signiE) 5{E^ - el) , (23) 

to a broad distribution of interacting particles, which is described by a double Lorentzian curve 
with half maximum width 2r: 

Tlo{E^+el + T^)-2E^-p + 2Em 

aHp, e) = - . (24) 

Here e^- denotes the onshell energy for relativistic fermions. Due to the non zero spectral width 
parameter F, the Fourier transform of (j24p into the time representation is exponentially damped. 
So the spectral width represents a finite lifetime of every excitation in the medium and can thus 
be associated with the microscopic scattering rate of the source particles. For F — > 0, one recovers 
the free spectral function A^p'^{p,E) from (j24p . 

For quarks in a QGP, the scattering width F has contributions both from electromagnetic and 
strong interaction processes 471 ]: 

F(r) = Fe^(r) + Fs(r) . (25) 

For the temperatures T, which occur in a quark-gluon plasma, the contribution from strong in- 
teraction processes is proportional to the temperature times some constant of order 1, i.e., it 



basically behaves as Fs(T) T [47|. So F may typically range from 0.1 to 1.0 GeV [4^, |49|]. The 



electromagnetic contribution behaves as T(.m{T) ~ aeT and can be neglected. 

Since we are only interested in the photon emission from a particle (not antiparticle) , we take 
into account only the fermion component of (j24p . which has its peak at positive energies: 



Upon transformation of (|2ip and (|22p into the time representation and evaluation of the trace in 
(dHI), we obtain for zH^j^(A;, t, u): 

«n<j,(A;,t,M) = — -{{PtJi{p-k)^+py{p-k)^-g^,y[p-{p-k)-rin?]] 



xe 



-i{ep~e^_^){t~u)^-2r\t-u\ ^27) 
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The superscript ' has been omitted for 



Here we have introduced po = £p and {p — k)o = £^_f^- 
convenience. If we now insert (j27p into ([8]), the volume factor V drops out when integrating the 
photon spectrum over d^x and we obtain 



dtdFk 



it) 



Pfi{p - k)u + Pu{p - k)^ - g^,y [p- {p-k)- m^] 



(2^)3 

X / du e-2r|*-"lcos(a;fc(t 

J —oo 



where tok is given by 



u)) 



k . 



(28) 



(29) 



It is helpful to rewrite (j28p in terms of A;|| and k±, which denote the components of the photon 
momentum parallel and perpendicular to the jet momentum respectively. We then carry out the 
integration over the azimuthal angle (pk and obtain 

kf 



k d^ 



k± dtdk\\dk_ 



vr 



1 



1 + 



p{p-h) 



m 



+ 2 




p{p — k\\) + m? 



+ 



2k^, p-ku 



k^ 



du e 



-2r|t-?i| 



cos {oJk{t — u)) 



(30) 



Expression ()30p allows for the investigation of the radiative behavior of a fermion jet in a non- 
equilibrated hot plasma. The expansion and cooling down of the plasma results in a time dependent 
scattering rate T = T(T{t)) = r(r). Due to the time integral in (j30p . the photon production rate 
at time t does not only depend on F at time t but also on F at earlier times. Therefore potential 
memory effects are covered by this expression. 

For completeness, we briefly discuss the applicability of the one loop approximation. It corre- 
sponds to the direct radiation from the different scattering points whereas interference contributions 
are described by so called ladder diagrams 50(]. So interference contributions have not been taken 
into account explicitly. Nevertheless, the ladder diagrams can be summed up to an effective one 
loop diagram when assuming that a nonrelativistic velocity of the source particle is reduced by 
a constant fraction at each collision [s^]. In this case, the summation procedure effectively leads 
to a reduced scattering rate F. Therefore, we may consider our one loop diagram as an effec- 
tive one where interference contributions and thus the LPM effect are already taken into account. 
Furthermore, the one loop approximation is formally the leading order contribution in the QED 
framework. 

One easily verifies that (j30p cures the infrared divergence of the standard pQFT-result [24. l25l]. 
This is because the problem of multiple scattering has been treated using full Green's functions 
with damping. The damping results from repeated scattering and leads to a suppression of the 
photon production rate for Uj^ ^ 2F. This range will from now on be referred to as LPM region. 
In the Bethe Heitler region, i.e., for uJk ^ F, our result reproduces (compare [50(] for the non 



relativistic case) the standard pQFT expression [24 . |25I | : 



d' 



n 



d^xdk 



1 

2^ 



2aePo^ 
k 



lij{k)(3i(3j 
{I -13 cos Of 



(31) 



Here (• 



denotes the integration over the solid angle in fc-space. We shall give an explanation 



of how the suppression mechanism for uj^ <^T comes about. The current of the source particle 
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decorrelates within the time interval ts = 1/{2T). As this is the timescale over which the source 
particle 'remembers' a scattering, radiation induced by this scattering is only possible within this 
time interval. In the LPM region, however, the formation time 

TF ^ 2TT/uJk 

of the photon is much larger than this decorrelation time. So the probability for the photon to 
be emitted within the decorrelation time interval is very small and the photon production rate is 
suppressed. This shows that the LPM effect is not purely an interference phenomenon. In fact, a 



perturbative treatment of interference contributions fails to cure the infrared divergence [51 1 



III. NUMERICAL INVESTIGATIONS AND RESULTS 

In order to investigate the radiative behavior of a fermion jet in a non equilibrated hot plasma, 
we simulate the evolving medium by varying the scattering rate T. The size of possible memory 
effects is most effectively extracted when changing F linearly over different time intervals At. Since 
the integrand of (j30p depends on two different time arguments we have to decide on which time 
argument r the scattering rate T shall depend. It is possible to introduce T = T{T{u)) = T{u) or 
r = r (r (^-^)) = r (^-^)), l-e., r can depend on the backward time u or the center of mass time 
The choice T = T[T{t)) = T{t) is also possible but it results in an entirely Markovian time 
evolution of (|30p . For our numerical investigations, we consider a p = 20 GeV fermion jet. 



^2 

time t 

FIG. 1: Linear change of F over At. 




As initial and final scattering rate, it is convenient to choose Fi = 1.0 GeV and F2 = 0.3 GeV, 
respectively. To find the timescales in which the photon rates adjust to changes in the medium and 
how these so called memory times depend on the momenta of the emitted photons, we consider the 
relative difference between the final (at t — > 00) and an intermediate photon rate which we define 
as follows: 

(fn{t ^ 00) - (fn{to + At) 



r(At) 



(32) 



d^n(t 00) — (Pn{to) 

Here to is the initial time at which F starts to change and to + At is the point of time where the 
change of F is complete. In the Markovian case, the photon production rate at to + At has already 
evolved to its final value so that r(At) = 0. In the case of full memory, the photon production rate 
at to + At still remains at its initial value and we have r(At) = 1. In general, it has a finite value 
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different from 1 and thus quantifies the memory effects for a given change duration At. Therefore, 
we are able to determine the memory times of the system by parametrizing r(At) and adjusting 
the parameters to the result. 

Figs. [5] and El show the relative difference r(At). For k± <C fcy we find an exponential decay 
r(At) ~ e""^*/^ (Fig. [2]) while for larger k± we have r{At) ~ (Ai/r)e~'^*/^ (Fig. E]). In both 
figures, r(At) is shown for F = T{u). 
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FIG. 2: r{At) for /cy 
GeV (LPM-region) 
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0.5 FIG. 3: r{At) for 
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So the following parametrization incorporates both cases 



r{At) 



T 



-At/T 



(33) 



The memory time is given by the parameter r, which describes how fast the system loses information 
about the initial state. Its values are shown in Figs. H] and [5] for fixed k\\ and varying k±. We 




FIG. 4: Memory times for F = F(m) FIG. 5: Memory times for F = F[(t + u)/2] 

find that the memory times become smaller at larger k±, a behavior that can be understood when 
regarding the characteristic timescale (r) that contributes to the memory integral in (j30p 



(t) ~ Re I ^—22 ^ i K = Zl f34^ 
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where ujk reads in terms of and k_\_: 



p^ + kl + k\- 2pk\\ - Jkl + kl 



(35) 



So the larger /c^, the shorter the contributing timescale to (j30p and hence the smaller the memory 
effects. In fact, if we compare Figs. 0113 with (|34|) we see that the memory times essentially behave 
as the characteristic timescale (r) which can be understood as follows. For kx_ <^ /cy, i.e., for 
LUk ^ 2r, the decorrelation time ts = 1/2T sets an upper bound for the memory time. That is 
why r only depends on Fi and F2 but not on k\\ in that range. One can indeed show [52| that the 
memory time behaves as r ~ I/F2 for k± 0. 

In the Bethe-Heitler regime {k± ^ the memory time behaves as t ^ tf/N where N is the 
number of oscillations in the integrand of ()30p during the time interval ts- So the expression that 
the memory time scales with is much smaller than the formation time Tp- Due to the oscillations 
only the fraction 1/N of the information from the timescale Tp remains. 

So the memory time r is not identical to the formation time Tp which one might expect. In 
particular, both quantities do not coincide with each other for k± <^ fcy ^ p. In this limit, Tp can 
be approximated as: 

k\\ lcs{k\i,k±) 
"^"fci= k^ ^''^ 

which is also known as Dokshitzer's approximation [5^. One could already infer this from the 
decorrelation time ts which sets an upper bound for the memory time whereas (j36p diverges for 
k± — > 0. For completeness, we shall briefly point out how (j36p comes about. In a reference frame 
where fcy = 0, the formation of a photon is given by Ti;' ~ 1/A;_|_. If we consider a Lorentz boost 
into a reference frame where /cy 7^ 0, the formation time obtains a dilation factor of approximately 
7eff(/c||, /c_l) ~ k\\/k± for k\\ » k±. 

We shall briefly explain why the memory times are only half as big for F = F (^-y^) as they are 
for F = F(n). At the time t, when the full change is reached, F = F (^-^) is equal to F = F(n) for 
twice the switching time interval At which can also be inferred from Fig. [6j So for F = F (^-2^), 



the photon rate after At has already developed as far as for F 
effects are only half as big. 



F(n) and At' = 2At. So memory 



r(u) 



r=r^,(u) 

r=r^,[(t+u)/2]=r2^,(u) 




FIG. 6: At the time t when the full change is reached, F = F (^^) for At corresponds to F = F(t 
for 2At. 



Even though both implementations render a causal time evolution, F = F (^-^) can nevertheless 
be considered as the more physical one. As the photon self energy in<^(A:, t,u) is the current- 
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current-correlator for the different points of time t and u, one intuitively expects that the change 
of r enters in^^,(fc, t, u) via both time arguments and not only via one of them (see Fig. [7|). Since 
the memory times are clearly below 1 fm/c for F = F (^-5^), our results already indicate that the 
adjustment process of the radiative behavior of the fermion jet to changes in the medium is almost 
instantaneous. 



(t + u)/2 




FIG. 7: The time dependence of F enters ill^j^{k,t,u) via both time arguments as this expression 
is the current-current-correlator for the points of time t and u. 



Now we turn to the total photon rate and the total radiation power, which are obtained by 
integrating (j28p over all photon modes. Since F does not depend on the photon momentum k, 
the total photon rate and the total radiation power are logarithmically and linearly divergent, 
respectively. This unphysical behavior is cured by introducing a cutoff at k = p, preventing the 
fermion jet from emitting photons with higher energy than its own: 

- f\^dk[dn,^^ (37a) 



dt Jq J dtd^k 

"^C* - (37b) 



dt Jq J dtd^k 

The memory times for these two quantities are calculated in the same way as for the individual 
photon modes and shown in Figs. [8] and [H Again the memory effects are only half as big for 
F = F (^^) as they are for F = T{u). The overall behavior corresponds to that of the photon 
modes in the LPM-region as these modes have significantly larger memory times than those in the 
BH-region. The memory times of the total radiation power are smaller than those of the total 
photon rate, because soft photons with larger memory times are suppressed by a factor of k. Again 
the adaption of both quantities to the scattering rate F is almost Markovian (r < 1 fm/c). 

Finally, we calculate the energy loss of the fermion jet during the time At. A comparison of 
our result to a quasistatic calculation, as shown in Figs. [10] and [TT\ reveals a difference of less 
than 5 percent for the characteristic expansion time {At = 4 fm/c) and scattering rates (Fi = 1.0 
GeV, F2 = 0.3 GeV) of a quark-gluon plasma. A significant aberration (> 10 percent) can only be 
observed for At = 1 fm/c. This shows again that the adjustment process is close to Markovian. 

In the case of At = 4 fm/c, we get an energy loss of —AE = 6.32 GeV which amounts to 
about one third of the fermion jet's onshell energy. Such a high energy loss is clearly unrealistic 
for photon emission even though the scattering rate is dominated by strong interaction processes. 
If we considered gluon emission instead, the electromagnetic coupling constant of ae ~ 1/147 in 
()30p would be replaced by the strong coupling constant of ag ~ 0.3 such that the energy loss would 
increase by a factor of roughly 40 and the fermion jet would be damped off immediately. 
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The main reason for such a high energy loss is the slow convergence of (I28p for large k. Without 
the cutoff in equations (|37ap and ()37bp the total radiation power and thus the energy loss would 
even be infinite. This behavior results from our ansatz for and 5^ which allows for transitions 
between arbitrarily offshell energy states which is not physical. In order to see this, we consider 
a modified ansatz where the fermion can only undergo transitions between energy states which do 



7 
6 
5 

> 4 
(3 

m 3 



+ r=r{u) 

X Markovian Case 



,X' 



2 
1 

ok 





,x 



,x 



10 
t [GeV"''] 



15 



20 



7 
6 
5 

¥ 4 
^ 3 



2 
1 

0* 




+ r=r[(t+u)/2)] 
X Markovian Case 



5 10 
t [GeV"''] 



15 



20 



FIG. 12: Energy Loss for At = 4 fm/c 



12 



not differ from the onshell energy by more than nT: 

S<{p + k,E) = 27ri^^6'-^\p + k-if)C{n)e{nT-\E-e^^f:\)AF{p + k,E) (38a) 

S^{p,E) = -2TTie{nT -\E -ej^\)AF{p + k) (38b) 

Here n can be any positive number and must not necessarily be an integer. The expression C(n) in 
Sp{p + k, E + k) is a normahzation factor which is given by C(n) = 7r/2 arctan(n). It is necessary 
to keep the normalization condition for a single spin- 1/2 fermion 

-^Vj ^Tr{^oS^ip,E)}=2. (39) 

The propagator Sp{p, E) is normalized to 2 because the spectral function only includes the particle 
and not the antiparticle component. One easily verifies that with the above prescription one has 
no more photon emission for cj^ > 2nT. Therefore the total radiation power and the total energy 
loss is tremendously reduced for small n which is shown in Table |l] for the parameters Fi = 1.0 
GeV, r2 = 0.3 GeV and At = 4 fm/c. The case n = oo corresponds to the original result. 



n 


dE/dt [GeV^] for Ti = 1.0 GeV 


dE/dt [GeV^] for T2 = 0.3 


-A^ [GeV] 


0.5 


0.017 


0.005 


0.22 


1 


0.052 


0.016 


0.70 


2 


0.120 


0.038 


1.60 


5 


0.239 


0.080 


3.25 


10 


0.327 


0.114 


4.53 


oo 


0.407 


0.199 


6.32 



TABLE I: Radiation power and energy loss for different values of n 



If we compare the energy loss —AE for different values of n and take into account that its value 
is scaled up by factor of ag/ae ~ 40 for the case of gluon emission, we see that our results become 
more realistic if n is chosen sufficiently small. For the case of n = 0.5, the energy loss increases to 
—AE = 8.8 GeV for gluon emission. This is almost 50 percent of the onshell energy of the fermion 
jet which from the phenomenological point of view seems to be quite reasonable. 

The energy loss in Table |T] has been calculated for the case of a Markovian time evolution. 
Within the original ansatz, the strong contribution of hard photon modes with small memory 
times is also the reason for the small overall memory effects. As the emission of hard photons is 
suppressed by restricting the number of available energy states, we investigate if this procedure 
also results in larger memory effects. For this purpose, we dynamically calculate the energy loss 
for n = 1 and n = 0.5 in Eqs. ([38]) and again compare our results to a quasi static calculation. 

We can infer from Figs. [13] and (TJ] that the memory effects are clearly larger than those for 
the original ansatz. We can now observe significant aberrations between the quasi static and the 
dynamic result for the characteristic expansion time of a quark-gluon plasma (At = 4 fm/c). They 
reach up to 15 percent for n=1.0 and increase even further (up to 25 percent) for n=0.5. So 
restricting the number of the available initial and final energy states does not only result in a 
tremendous reduction of the energy loss —AE but memory effects also play a sizable role as the 
contribution of hard photon modes with small memory time is suppressed. 
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But we have to keep in mind that the size of the memory effects naturally depends on n and 
that we are principally allowed to choose this parameter freely as the energy loss —AE remains 
finite for any finite value of it. This in turn results in smaller memory effects for larger values of n 
which is depicted in Figs. [T5]|17i So getting more reliable results about their relevance for radiative 
jet energy loss probably needs a full microscopic and quantum field theoretical approach to this 
phenomenon. 



IV. SUMMARY AND OUTLOOK 



In this work, we have investigated the electromagnetic radiative energy loss of a fermion jet 
which moves through the evolving quark-gluon plasma created during heavy-ion collisions at high 
energies. We focused on whether the radiative behavior of the jet adjusts to the evolving medium 
instantaneously or if existing memory effects are of importance. For this purpose we have applied 
the real time formalism of finite temperature field theory, which allows for a full quantum mechan- 
ical treatment of such a non-equilibrium phenomenon. Using an effective in-medium propagator 
for the fermion-jet, we determined the time scales in which the radiative behavior of the jet adjusts 
to changes in the medium, and also calculated the energy loss of a fermion jet over the time in 
which the plasma expands and cools down. As a first step, we have restricted ourselves to photon 
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emission. 

For the characteristic parameters of a quark-gluon plasma (effective initial and final scattering 
rates of the jet quark: Fi = 1.0 GeV, and F2 = 0.3 GeV and a life time of approximately At = 4 
fm/c), the memory times for the evolving total photon rate and total radiation power have been 
found to be in physical terms clearly below 1 fm/c. Furthermore, a comparison of our result for 
the total energy loss with a quasistatic calculation has revealed a difference of less than 5 percent. 
This supports the validity of quasi-equilibrium calculations when calculating photon production 
from a QGP and indicates that a full dynamic treatment is not crucial. The main reason for 
that are the large scattering rates in a QGP leading to a fast decorrelation of the source particle 
defining the major time scale for memory effects. If one considers dilepton production from the 
subsequent hadronic phase where the scattering rates are smaller 35|, ISg] , memory effects play a 
more important role. 

But for completeness, it is important to point that the small memory effects also result from our 
specific ansatz for the fermion propagators which allows for transitions between arbitrarily offshell 
energy states. This results in a strong contribution of hard photon modes with small memory 
times leading to the small overall memory effects mentioned above. As our ansatz requires the 
introduction of a cutoff in order to keep the total photon rate, the total radiation power and the 
total energy loss finite, it is strictly speaking not physical. For that reason, we have also considered 
a modified ansatz where the number of available energy states was restricted to a finite interval 
around the onshell energy which reduces the emission of hard photons. Besides a strong reduction 
of the totally emitted energy —AE, this also leads to significantly larger memory effects. The 
magnitude of these effects does, of course, depend on the chosen width of the energy interval. 
Thus, even though we were able to determine the role of memory effects for all photon modes 
separately, finding a final answer to the question whether they are relevant for the total radiative 
jet energy loss might indeed require a full microscopic and quantum mechanical treatment. 

In this work, we have restricted our considerations to QED processes. The extension of our 
investigations to the non-Abelian QCD will bring complications. Gluons can interact with each 
other and additional self energy diagrams become relevant. Another important problem in the 
QCD framework will be that the formation time of the emitted gluon depends on the momentum 
transferred to the gluon during the rescattering in the medium. This phenomenon has already 
been investigated in [28 1. 
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